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ABSTRACT 

Web-stiffened sandwich cylindrical shells consisting of 
two concentric cylindrical shells joined and stiffened by 
annular rings (webs) were considered in the investigation. 
The purpose of this investigation was to study the stress 
distribution and buckling modes, and to determine the scant- 
lings and configuration of optimum proportions. 

The scantlings and configuration of optimum proportions 
was defined as that with the lowest weight to displacement 
ratio at a given loading of hydrostatic pressure. 

The scantlings (outer shell thickness, inner shell 
thickness, web thickness, web depth (separation of shells) 
and configuration (web spacing) were systematically and 
independently varied and the effect upon failure mode pres- 
sures, and the weight t<3 displacement ratio was determined. 
The results were plotted graphically. 

The form of stress and buckling analysis results follow 
closely those of a ring-stiffened cylinder. Geometric sim- 
ilitude was determined to exist. The scantlings and configu 
ration of optimum proportions were found to be a function of 
the loadings. 

Web-core sandwich configurations possess structural 
efficiencies on the order of 20% higher than those of con- 
ventional ring-stiffened construction for loadings on the 
order of 1000 psi (2225 feet) and these diminish with depth 
to about 5% higher for loadings on the order of 4000 psi 
(8900 feet) . 

The real advantage of the web-stiffened sandwich con- 
struction lies in the use of thinner plating with all its 
superior metallurgical and fabrication qualities. 

Thesis Supervisor: Professor J. H. Evans 

Title: Professor of Naval Architecture 
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NOTATION 



b 


Web thickness 


c 


‘W - 1 


C 


Web constant 


D 


Bending, rigidity, Et 3 /12(l-v 2 ) 




Mean diameter of shells 


a 


RRO-RRI 


E 


Young 1 s modulus 


E 

s 


Secant modulus 


e t 


Tangent modulus 


f 


Stress ratio, a /a, 

X (p 


*o' *i 


Edge coefficients for outside and 
inside shells 


^oo' ^oi' ®io' 


a.. Lame deflection coefficients for the 


(M,N,P,Q) 


annular webs 


h 


R-R. 

O 1 


k. 

i 


Stress intensity per pound of pressure 
(K x = V? 1 


L 


Unsupported length of cylinder 


l b 


Length between bulkheads 


m, n, k 


Integers 


N x' N * 


Longitudinal and circumferential 
forces per unit length 


X, <p, r 


Longitudinal, circumferential, and 
radial coordinates 



B . , F. , G. 

10 10 10 


Hyperbolic and trigonometric stress 
parameters for outside shell 


B il' F il' G il 


Hyperbolic and trigonometric stress 
parameters for inside shell 


M x 


Longitudinal moment per unit length 


p 


External hydrostatic pressure 


g , d , f 
o o 


Force coefficients for outside shell 


H 

H 


Force coefficients for inside shell 


H , H. (HHO.HHI) 


Total radial forces on the outer and 
inner surface of webs per unit length 


V W v i> 


Longitudinal forces on outside and 
inside shells per unit length 


Qq r Qi ( Ho / i ) 


Radial forces on outside and inside 
shells per unit length 


( ) 


Same quantity used in general theory 

and computer programs 

% 




Radial force per unit length 


R 


Radius 


R , R. 
o' 1 


Mean radius of outside and inside 
shells 


RRO , RRI 


Extreme radius of outer fiber of out- 
side and inside shells 


RR 


RRI /RRO 


S 


Mean distance between webs 


t 


Thickness of shell 


T , TO 


Outside shell thickness 


T i# TI 


Inside shell thickness 
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U, V, w 


Longitudinal circumferential and 
radial deflections 


W , W ■ 
os ' is 


Radial deflection of outside and 
inside shell 


W , W • 
ow ' 1W 


Radial deflection of the outer and 
inner surfaces of the web 


Z 


Radial bending rigidity, E/(l-v 2 ) 


V V £ r 


Longitudinal, circumferential, and 
radial strains 


°x' V 0 r 


Longitudinal, circumferential, and 
radial stresses 


CD 

0 

CD 


Shell flexibility parameters for out- 
side and inside shells, 

[ 3 ( 1-v 2 ) ] VC [L/(Rt) 1 / 2 ] 




3ending stress coefficient, 

[ (1-v 2 ) /3] V 2 


\) 


Pbisson's ratio for isotropic 
materials 


a . 

i 


Stress intensity, Huber-Hencky , 
Von Mises , ( - o-^a 2 ) 1 /4] 


a 

Y 


Yield stress 


u ijKL ^°iKL^ 

i = x, <£>, r 


Stress location 
direction 


j = F,M 


shell (web-shell juncture, midbay) 


K = 0,I,W 


outside shell, inside shell, web 


L = 0,1, M 


fiber (outer, inner, membrane) 


(jiKL) i. 


. e. , MXOI = midbay, longitudinal, outside 
shell, inner fiber 



P 



Yield failure at shell midbay 



PF Yield failure at web-shell juncture 

PWS Yield failure at web outer fiber 

PPC1 Axisymmetric shell buckling 

PPC2 Asymmetric shell buckling 

PPC3 General instability 

PPC4 Web, ring buckling 

Web, edge compression buckling 



PPC5 



INTRODUCTION 



One of the more promising structural concepts for the 
design and fabrication of cylindrical pressure hull struc- 
tures is the sandwich concept. Structural engineers in the 
aircraft industry have long recognized and taken advantage 
of the favorable strength-weight characteristics of 
sanawicn-type construction. In studying the literature, 
however, it nas been found that the loading conditions 
encountered in tnese applications nave dictated sandwich 
structural arrangements that would be of no direct use in 
the design of pressure hulls for submersibles . 

The demands of hydrostatic pressure loading are such 
that in order to exploit tne sandwich concept for pressure 
hull construction a core comprised of elements which are 
compression resistant as well as shear resistant is sought. 
The major difference is that membrane loads are predominant 
in hydrospace applications , wnereas in aerospace, the bend- 
ing and shear type loads are of prime concern. 

Experimental programs and analytical studies nave been 
going on concurrently in order to develop rational formulas 
based on thin-shell theory for predicting the static struc- 
tural response of these types of structures. In reference 
10 Pulos presents an analysis of the axisymmetric elastic 
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deformations and stresses in a web-stiffened sandwich 
cylinder under hydrostatic pressure. Raetz presents a 
similar analysis for the toroidal tube-stiffened sandwich 
cylinder . ^ In references 12 and 14 Nott presents a sim- 
plified stress ana strain analysis. 

Nott ' s equations were used and programmed for the 
stress analysis of the web-stiffened structure. 

In reference 22 the axisymmetric and asymmetric buck- 
ling equations are developed. In references 15 and 23 an 
expression for the general instability is determined. 

The above equations are utilized, with others derived 
for web stability, to study the parametric nature of web- 
core sandwich structure and to determine optimum propor- 
tions . 

* 

Optimum design must be based on rational considera- 
tions of inelastic behavior, however, and not on the "one 
horse shay" concept based on elastic considerations of in- 
stability. The "ignorance factors" can then be representa- 
tive of the variability introduced by certain intangibles 
which are not easily considered in a theory. 

Some of the intangibles which influence static strength 
and complicate the problem so that appropriate design for- 
mulas cannot be derived on purely theoretical grounds are 
those inherent in the fabrication process itself: 
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(1) Initial stress. 

(2) Residual stress. 

(3) Imperfect circularity. 

(4) Non-isotropic or non-homogeneity of material. 

(5) Actual boundary conditions at stiffening rings. 

Therefore, in view of these complications, the devel- 
opment of satisfactory design criteria must first be pre- 
dicted on rigorous mathematical theory with its concomitant 
idealizations, and then, empirical factors derived from 
test data can be introduced to "adjust" the theories to 
take account of these many variables which could not or 
were not considered. 

Thus an inelastic analysis of all instability modes is 



used. 
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GENERAL THEORY 

1 . Stress Analysis 

A theoretical analysis of the axisymmetric elastic 
deformations and stresses in a web-stiffened sandwich 
cylindrical shell structure under external hydrostatic 
pressure is presented in Appendix A. The solution is based 
on the use of edge coefficients for plate and shell ele- 
ments of finite length, and includes the computation of 
the edge forces and moments arising at the common junc- 
tures of these elements. 

Equations are given for computing numerically the 
longitudinal and circumferential stresses in the two 

coaxial cylindrical shells and the radial and tangential 

* 

stresses in the web stiffeners between the two shells. 

To evaluate the structural strength of cylindrical 
sandwich shells, the locations and magnitudes of maximum 
stresses must be determined. Rigorous analyses for stress 
distributions in sandwich shells loaded under external 
hydrostatic pressure were carried out by Pulos^ and 
Raetz J " L . These analyses illustrate that the maximum 
stresses occur in the inner and outer shells at locations 
next to the annular webs and midway between the two webs 
that separate the shells. Therefore, the stresses in the 
inside and outside shells at these two locations are of 
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interest for the purpose of structural design and 
evaluation. 

The procedure and nomenclature for programming the 

12 14 

stress analysis is that of Nott. ' 
a. NOMENCLATURE 
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b. ASSUMPTIONS 

From symmetry considerations it is seen that the 
edges of the web stiffener, circular annulus, do not under- 
go any rotation. This stems from the fact that a horizontal 
tangent or zero-slope condition is assumed to exist at the 
junctions of the webs with the two cylindrical shells. 

This assumption implies that the edge moments on each shell 
at the shell-web juncture balance each other, so that there 
are no net moments to be resisted by the web. Further, it 
is assumed that the web elements do not take any axial 
force due to the axial pressure, but that this is all re- 
sisted by the cylindrical shells. 

Also inherent in the analysis (see Appendix A) is the 
neglect of the beam-column effect due to the axial portion 
of the hydrostatic pressure. 
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D 



d 4 w 
dx* 



+ N 



d 2 w 



x dx 2 



Et „ 

+ — T W = 
R 2 



- * N 
R 



x 



neglecting the beam-column effect becomes 



D 



d **w 
dx 1 * 




P 




x 



Thus the analysis of the web stiffeners reduces to 
that of a circular annulus subjected to axisymmetric in- 
plane radial forces on both its inner and outer boundaries. 

On the basis of these assumptions, it is only necessary 
to derive edge coefficients for an annulus undergoing radial 
deflection. 

The functions A, N, P, and Q are the Lame deflection 
coefficients for the annular webs and can be expressed 



M = 


Rw 
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(1 + 




d 

2R W 


2v ) 


N = 


Rw 

b 


U + 2Rw } 




d 

2R W 


2) 


P = 


*w 

b 


( 1 2rJ 




d 

2Rw 


+ 2) 


Q = 


^w 

b 


U 2Kw’ 




d 

2R^ + 


2v) 



and the functions g Q , g^ are the edge coefficients from the 
outside and inside shells: 
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where 
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2R ^ 



T o L 
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10 



2R? 
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T ± L *li 



10 
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o o 
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0 O 

sinh 0 . + sin 0 . 
i i 

_cosh 0. - cos 0. 

1 l 
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✓3 (l-v2) 
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o o 




✓ 3'(l-v2) 



/s-f- 



Thus, from Appendix A, the axial and transverse forces 
(V 0 , , H o and H^) , which act at the intersections of the 

annular webs and cylinders can be determined. 
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And from these forces the maximum stresses in the 
shells can be determined. 
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C. MAXIMUM STRESSES 
OUTER SHELL 
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and likewise for B^, 32^/ . F^, ky substituting the 

inside I for the outside O. 



d. WEB STRESSES 



From Appendix A the maximum stress in the web 
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rwo 



[ (-HHI (RR) 2 +HHO) - (-HHI ( RR) +HHO ( RR) ) (RR) ] 
b/2 (l-(RR) 2 ) 
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4>wo 



[ (-HHI (RR) 2 +HHO) + (-HHI (RR) +HHO (RR) ) (RR) ] 
b/2 (l-(RR) 2 ) 



where 

HHO = H Q + b/2 
HHI = H. 

l 



RR = RRI/RRO 
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2 . Buckling Analysis and Failure Modes 

The axisymmetric and asymmetric (lobar) buckling 

; 4 g 

equations developed by Lunchick and Reynolds for ring- 

stiffened cylinders are modified in Appendix 3 for applica- 
tion to web-stiffened sandwich cylinders. They are eval- 
uated for both the inner and outer shells. Reference 22 
indicates good agreement between theoretical and experi- 
mental collapse pressures. 

The general instability equation used has found con- 

15 23 

siderable experimental verification. ' 

The two modes of web buckling are first, a variation 
of the Foppl or Levy ring buckling (effectively TOKUGAWA) 
and second, the buckling cf an annulus in biaxial edge 
compression (see Appendix B) . 

a. FAILURE MODES 
1. STRENGTH 

Using the Huber-Hencky-Von Mises criterion 



c . = /a i + o 2 - a , a 
l ip x cj> x 



a • 



i - P ^ K | + k£ - K^K 



<J)^X 




+ K £ - K^K X 



evaluated at the web-shell interface (PF) , midbay (P) , and 
outer fiber of web (PWS) . 
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2. AXISYMMETRIC BUCKLING 
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3. ASYMMETRIC BUCKLING 
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4. GENERAL INSTABILITY 



PPC3 = SE 



T *V 



n 2 -1+ A - 



2 1 



(n 2 +A 2 ) < 1-3/41 l-= 



1 



J T 



( n 2 +A 2 ) 2 



+ 



(n 2 -l) E t I eff 
r 3 l 



22 



where 
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5. WEB BUCKLING 



a« Ring Buckling 
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b . Edge Compression 
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AXI SYMMETRY BUCKLING ( BETWEEN WEBS) 




A SYMMETRIC. (LOBAR) SUCKLING (BETWEEN WEBS) 




GENERAL INST ABU. ITT ( CONCURRENT FAllU RE 
OF SHELLS AND WEBS) 

FAILURE nODES (SHELL) 




we b ( ring) suckling 




WES ( com prsssion) suckling 



FAILURE Njopgs ( W E 6) 
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PROCEDURE 



The buckling pressure is a function of the cylinder 
geometry and the secant and tangent moduli as determined 
from a stress-strain intensity diagram for the shell 
material . 

Before the buckling equations can be used, E g and E^ 
must be related to the applied pressure. The secant and 
tangent moduli are defined 




E t = da i /de i 



and are shown graphically: 
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Hence E and E m are readi.lv determined from the stress- 
s T 

strain curve (see Appendix C) . Having determined E /E and 
E t /E as functions of a ■ , they must be applied to the hydro- 
statically loaded cylindrical shell. By expressing hydro- 
static pressure .in terms of the stress intensity, a rela- 
tionship between E , E„, and pressure will be established. 

o i. 




G\ 

(1) P c = f (geometry, E g , E^) 
a buckling equation 

(2) a i = SoJ + - a 1 a 2 = p/K* + K* - 

a stress intensity (Huber-Hencky-Von Mises) 
equation 

The intersection of these two curves will predict fail- 
ure by inelastic buckling if below the yield strength. 



Thus the general solution technique is to (1) solve 
for the state of stress, and consequently a relationship 
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between and p, and then (b) solve the buckling modes as 
a function of o.(e /E^) . The intersection of the curves 
being the critical pressure (see lb) . 

1 . Solution Technique , 
a. STRESS ANALYSIS 
INPUT DATA: 

Rq / R j / / T j / b / L/ V/ p 

Let p = 1.0 psi 

(1) Formulate 6 , 8^ 

(2) Formulate F 1Q/ F^, F 2Q/ F 21 ... 

(3) Formulate g Q , g^ 

(4) Formulate , d 

(5) Formulate M, N, P, Q 

(6) Formulate V q7 

(7) Solve simultaneous equations for H q/ H^ 

(8) Formulate HHO, HHI , R R 

(9) Formulate a uQ , a ui 

(10) Formulate B 1Q , B-^, B 2C)/ B 21 ... 



(11) Solve for the states of stress 



b. FAILURE MODES 
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(1) INELASTIC MATERIALS 





(2) ELASTIC - PERFECTLY PLASTIC MATERIAL 





c. OPTIMUM PROPORTIONS 



A simple optimization program was used to deter- 
mine optimum proportions. As indicated in the flow 
diagram, a comparison is made between the design col- 
lapse pressure (PCO) and all the failure mode pressures 
(P^) . The failure mode pressures consist of: 

P E yield at midbay 

PF E yield at web-shell junction 

PWS E yield in web 

PPC1 E axisymmetric buckling, each shell 

PPC2 e asymmetric buckling, each shell 
PPC3 e general instability 
PPC4 e web instability (ring) 

PPC5 E web buckling (edge compression) 
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RESULTS 



1. Geometric similitude was determined to exist. The 
stress analysis and failure mode programs output are in 
Appendix D for three geometrically similar structures. 

Graph 1 shows the failure modes for these three geometric- 
ally similar structures made from HY-160 steel. The effect 
of the inelastic range of this material on the buckling 
modes is indicated. The example has underdeveloped webs 
(PWS, PPC4). 

2. Maximum stress states were found to be located as 
shown below. 



3. The results are plotted in graphical form. The results 
are divided basically to show (a) the failure modes of the 
shells, (b) the failure modes of the webs, and (c) optimum 
proportions and weight/displacement results. The presenta- 
tion of the results in this form was (1) to facilitate 
comparison with ring-stiffened cylinders, and (2) to 





MU 



show the relationship of the various failure modes. 

The form of the results are very similar to those of ring- 
stiffened cylinders. 

4. Graphs 2 through 7 are sample results of tne stress 
and buckling analysis for the shell portion of the web- 
stiffened sandwich cylinder using HY-30 steel. 

5. Graphs 8 through 11 are sample results for the web 
portion non-dimensionalized by P 

where 

D = 2 < T o + Ti> °y 
ND % 

The material was HY-80 steel. 

6. Graphs 12 and 13 show the effect of yield strength and 
the modulus of elasticity on shell failure modes respect- 
ively . 

7. Graphs 14 and 15 and Tables 1 and 2 are concerned with 
the optimization results. The optimum proportions for ring- 
stiffened cylinders shown in Graph 15 were based upon using 
92, 92A criteria for shell yield (outer fibers reaching 
yield stress) . 

8. For a given hydrostatic loading and hull diameter the 
shell thicknesses for the web-stiffened sandwich are approx- 
imately half that of a ring-stiffened cylinder. 



“Nl ’<JDS ^J3d SClNOOd d 3dOSS3>Jd SSdMIdOj 



GRRPH 1 




M2 128 l<H 

S T f\TE of stress (.G’O 



I Tt. 



140 



COLLAPSE pressure 



50 . 




.03 .OS .10 .2 0 .30 

UNSUPPORTED LENGTH -r DIAMETER (^/ dm ) 



.02 



.50 



COLLAPSE PRESSURE 'P POUNDS PER SQ. IN. 



51 . 




.02 .03 . oS .10 .20 .30 .SO 



UNSUPPORTED LENGTH -s- D1AMETER(%„) 



COLLAPSE PRESSURE P POUNDS PER SQ. IN. 



52 . 




o z .03 .OS -lO .zo .30 .50 

UNSUPPORTED LENGTH -7- Dl AMET ER (^DtO 



COLLAPSE PRESSURE P POUNDS PtR SQ. IN. 



53 . 




.03 .OS’ .10 .20 .30 .SO 

UNSUPPORTED LENGTH -r DIAMETER 



.02 



COLLAPSE PRESSURE P POUNDS PER SQ. IN. 



54 . 




i i t i 1 x I 1 i S 

.02. .03 .OS .10 .'ZO .10 .SO 



UNSUPPORTED LENGTH-.- Dl AH E TE R ( L y^ ^ 



COLlfcPSe PRESSURE P POUNDS PER SQ. IN. 



55 . 





56 . 



GRAPH 3 




WEB DEPTH -t SHELL AVERAGE THICKNESS , h/t s 



57 . 



QRPPH 3 




o 

in 



o 

6 



o 

in 






WEB DEPTH -T- SHELE AVERAGE THICKNESS , h/ t 



53 



GRAPH lo 




WEB THICKNESS-*-SHEU AVERAGE THICKNESS b/t s 



collapse pressure p pounds per sq. in. 



59 . 




.02 .03 .os .IO .20 .30 .50 

u N SUPPORTED LENGTH ~ Dl AMETER(V d^) 




UNSUPPORTED V. EtsiGTH -r Di AINETE R (%*) 



C OLL /\PSt PRESSURE P POUNDS PbR $Q. IN. 



61 . 




.03 .os .10 .20 .30 

UNSUPPORTED LENGTH -r DIAMETER y 



. 02 . 



.50 



qRRPH 




0 

1 



o 

o 

o 

00 



o 

o 

o 

r- 



o 

o 

o 

o 



o 

0 

0 

ta 



o 



o 

o 

0 

*0 




i 0 

! O 

_ o 

cst 



<3 S © 



o 



co 

0 



vS» <s/ 

0 0 C 



WDVO IldO 

HJ.c£3a 03^ lUN0iSN3U.ua -Mo N 



62 . 



2 

O' 

(/) 

or 

UJ 

a 

o 

D 

Z 

3 

O 

CL 



Cu 



COLLAPSE PRESSURE 



collapse pressure p Pounds Pt r Sq. in. 



GRAPH IS 



63 . 




64 



Table 1 

OPTIMUM SCANTLINGS 



P coll 


t/R 7 
W 


h/s 


500 


. 053 


1.43 


1000 


.095 


1.56 


2000 


.113 


2.14 


5000 


. 148 


2.16 



P coll 


h/R 
7 w 


h/s 


500 


.077 


2.50 


1000 


. 077 


1.87 


2000 


.077 


1.73 


5000 


.113 


2.14 



HY-30 



h/t s 


b/t 

s 


t /D 
o' 


17.65 


.583 


.0019 


17.39 


. 696 


.003 


12.00 


. 600 


.005 


6.15 


. 538 


.014 



h/t 

s 


HY-150 

b/t t /D 

' s o 7 


39.02 


.976 


.0011 


26.67 


.933 


.0016 


14.29 


.893 


.0028 


9.23 


.462 


.008 



t i /D 


L/D 


W/D 


.0015 


.020 


.137 


.00275 


.030 


. 239 


.005 


.025 


.436 


.014 


.030 


1.061 



V D 


L/D 


W/D 


.0010 


.015 


.125 


.0014 


.020 


.155 


.0028 


.020 


.258 


.005 


.025 


.522 



This value seemed high and thus a scatter check was con- 
ducted about this point. Results of this check are in 
Table 2. 
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Table 2 

OPTIMUM SCATTER ABOUT 









P coll 


= 500 


ps i 








P coll 


Mode 


VR* 


V b 


t i /D 


L/D 


b/t s 


W/D 


W/D* 


508 


PF 

o 


.058 


■ . 0014 


.0014 


.025 


.786 


.124 


.122 


554 


PWS 


.058 


' .0014 


.0014 


.020 


.643 


.128 


.115 


524 


PWS 


.058 


.0015 


.0014 


.025 


.759 


.127 


.121 


581 


PWS 


.058 


. 0016 


.0014 


.020 


.600 


.134 


.115 


647 


PPC4 


.058 


.0018 


.0014 


.020 


. 625 


.145 


.112 


677 


PPC4 


.058 


.0018 


.0015 


.020 


.667 


.153 


.113 



* sno 

Normalized W/D by ) . 

coll 



